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We propose a general framework that unifies the point of view of counting statistics of transmit- 
ted (fermionic) charges as it is commonly used in the quantum transport community to the point 
of view of counting statics of phonons (bosons). As a particular example, we study on the same 
footing the counting statistics of electrons transferred through a molecular junction and the corre- 
sponding population dynamics of the associated molecular vibrational mode. In the tunnel limit, 
non-perturbative results in the electron-phonon interaction are derived that unify complementary 
approaches based on rate equations or on the use of non-equilibrium Green functions. 
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I. INTRODUCTION 



The concept of counting statistics of transmitted 
charges as it was first defined and computed for non- 
interacting tunnel junctions by Levitov et al [l|, |^ is 
based on an analogy between coherent electronic trans- 
port and quantum optics [3| . In both fields of research, 
the underlying electronic or photonic quantum field un- 
dergoes quantum fiuctuations encoded in a fundamental 
quantity, namely the probability distribution P{q) that 
q charges are transferred through the junction or the 
stationary probability distribution P{n) of detecting n 
photons in the electromagnetic field, in a given amount 
of time. 

This concept has proved to be a powerful one [J| 
by providing deep connections between the fields of 
molecular electronics [5| and quantum noise [6|. On a 
more fundamental point of view, the concept of counting 
statistics of transmitted charges drew the attention of 
the condensed-matter community to the idea already 
present in quantum optics, that the quantum nature of 
transport mechanisms as well as the peculiar effects re- 
lated to quantum mechanics are more generally encoded 
into correlation functions of the fields rather than into 
the mean value of a given observable. 
Recently, ideas of measuring phonon shot noise have 
been reported in the literature [7|, but interestingly, 
there is to our knowledge no work that considers both 
kinds of descriptions on the same footing. We thus 
assign the following goals to the present article. It is 
first to provide such a general framework, that unifies 
the point of view of counting statistics of transmitted 
(fermionic) charges as it is currently accepted in the field 
of quantum transport, to the one of counting statistics 
of phonons (bosons) inherited from the quantum optics 
community. Then, it is to illustrate the fertility of 
such a framework by studying the concrete example of 
a quantum transport problem through an interacting 



nanoscale device. The corresponding issue is indeed 
relevant to the quantum transport community, where we 
can distinguish between two main approaches to derive 
the associated transport properties. 
On one side, the study of single electron transistors 
has motivated theoretical investigations of sequential 
tunnelling of charges through local devices, in presence 
of interactions with the electromagnetic environment [8| 
or with internal degrees of freedom (electron-electron 
or electron-phonon interactions). The corresponding 
theoretical framework based on rate equations with 
transition rates derived from the Fermi golden rule 
has been successfully applied in more recent works 
related to the problem of Franck-Condon blockade 
[9| and observed experimentally in suspended carbon 
nanotubes [10|. In presence of electron-phonon (e-ph) 
interactions, the obtained current versus voltage curves 



llj l exhibit characteristic patterns (satellite peaks at 



voltages multiple of twice the phonon energy) due to 
the activation of tunnelling assisted by phonon emission, 
as well as a reduction of the low-bias conductance 
(Franck-Condon blockade). The underlying transport 
mechanism was shown to induce giant Fano factors for 
the fluctuations of the transmitted charges (current 
noise), when the vibrational po pul ation was driven far 
into the non equilibrium regime |12| . The corresponding 
counting statistics of transmitted charges was derived 
for molecular junctions in Ref. llSJ, |lj| and for nanoelec- 
tromechanical systems in Ref. [l5[ . 

On the other side, the experimental works related to 
coherent transport through atomic chains [l6| and 
molecular junctions [l7| has driven a tremendous 
amount of theoretical work, based on the extensive use 
of Keldysh non-equilibrium-Green-functions formalism 
18 1 . In both approaches based on model Hamiltonians 



19 -21 1 or on more sophisticated ab-initio calculations 
22-|2J|, the conductance characteristics exhibit a jump 



at the inelastic threshold (voltages corresponding to 



the phonon energy) associated to the activation of 
tunneUing with emission of a phonon, the sign of which 
is determined by the transmission coefficient of the 
junction. Such feature is of great experimental interest 
in order to perform an inelastic spectroscopy of the 
device [20[ (spectroscopy of the vibrational modes, 
measurement of the e-ph coupling strength from the 
height of the jump of conductance). More recently, some 
theoretical works focused on inelastic signatures of e-ph 
interactions on noise characteristics 25|, |26| and on the 
full counting statistics of the transmitted charges in the 
coherent regime |27H29f . The appearance of jumps in 
the derivative of the current noise versus voltage curves 
was shown to result from a competition mechanism 
between elastic and inelastic processes of tunnelling, and 
the feedback of the phonon dynamics to be an important 
effect in order to correctly capture the behaviour of 
higher order cumulants in the regime of large voltages 

The organisation of this article is then the following. In 
the first part, we provide a general framework that uni- 
fies the point of view of counting statistics of transmitted 
(fermionic) charges in nanodevices to the counting statis- 
tics of phonons (bosons). As a particular case of our for- 
malism, we then consider in detail the particular model 
of a nanoscale junction as a single molecular level cou- 
pled to a local vibrational (phonon) mode and to perfect 
reservoirs with different chemical potentials (non equilib- 
rium situation). We show that both approaches based on 
rate equations in the sequential tunnelling regime or on 
non-equilibrium-Green-functions in the coherent regime 
may be described in a unified way following the scope of 
our formalism. We interpret the transport properties in 
such a device as the result of a self-consistent mechanism 
between charge transfer through the junction and the dy- 
namics of population of the local phonon mode. As an 
illustration, the full counting statistics (FCS) of trans- 
ferred electrons and phonons are defined and computed 
exactly in the tunnel limit. 



II. GENERAL FRAMEWORK 



H^^ = eVl2 




FIG. 1: Representation of the nanoscale junction. 



with a local phonon (vibrational) mode of energy loq. 
The electron-phonon interaction is proportional to the 
electronic density operator n^ ^ S d and to the vibronic 
position operator (f) — a + a^ . It is characterised by 
the e-ph coupling strength A. The second term in EqU] 
is the Hamiltonian of non interacting left (right) leads 
Hl{R) both maintained under a symmetric potential 
drop fiL = —fJ-B = V/2. In such a system, charge might 
be transferred from the electrodes to the molecule by the 
tunnel Hamiltonian Ht = (^Ld^^ + iiid^Jj) rf + he that 
couples the dot level to the X ^ L,R lead through the 
hopping term txd- The typical energy scale correspond- 
ing to this charge transfer mechanism is given by the 
X = L,R tunnelling rate Tx = \txd\'^/W, where in the 
wide band approximation, the bandwidth W = 1/ttpq is 
inversely proportional to the flat density of states of the 
leads pq. 

Although simple, the model described by Eq[T] exhibits 
a non trivial phase diagram characterised by two 
competing energy scales [l9|, namely F = F/, + F/j the 
total coupling strength to the leads and A^/wq the e-ph 
coupling strength. In this article, we concentrate on the 
regime of weak tunnelling for which F ^ A^/wq is the 
smallest energy scale of the problem. In this regime, the 
e-ph interaction strongly normalises electronic degrees 
of freedom and a non-perturbative approach in the c-ph 
coupling strength is needed. 



A. The model 



In the following, we consider the simple model of a 
molecular junction (see Fig[T]) described by the Hamilto- 
nian (in units e = h = k = 1) 



H — Hmol + 2_^ ^X + ^T 



(1) 



X=L,R 



In EqlU the molecular system is encoded into the 
Hamiltonian Hmoi = ('^d + A0)nc; + ujQa'^a that describes 
a single molecular (dot) level of energy LOd interacting 



We further perform a unitary Lang-Firsov (polaron) 
transformation [31] U — gandia-a ) -^^ order to explic- 
itly eliminate the e-ph interaction term from Hmoi by 
shifting the vibronic position operator (f> (g = A/wq is 
the dimensionless e-ph coupling strength) . The obtained 
dual representation is then more adapted for perturba- 
tion calculations in leading orders of F, and the trans- 
formed Hamiltonian H = UHW now reads 



H = H, 



mol 



X=L,R 



Hx + Ht 



(2) 



In Eql2l the Hamiltonian describing the molecular system 
Hmoi = Wdnd + woa^a does not couple anymore electronic 
and vibrational degrees of freedom, and the dot position 
is renormalised by the polaronic shift oj^ = uJd ~ g^^Q- 
Unfortunately, the price to pay in the Lang-Firsov trans- 
formation is that the transformed tunnelling operator 



Ht = {tLd^l 



tRdm^j^}dX 



he acquires a phase oper- 



ator X = e"*-" "" '' that controls the injection of charges 
in the system. 

B. Counting statistics of electrons and phonons 



namely ti(^)d ^ tL(R)def'^^*^'^ and a ^ ae-^«(*)/2. The 
mean value is taken with respect to the unperturbed 



Hamiltonian Hq = H — H' 



T;x{t),i{t), 



and the electronic 



counting field x(0 is equal to -l-(— )x on the positive 
(negative) branch of the Keldysh contour C. It stands 
for a non destructive measurement of the charge trans- 
ferred from any electrode to the dot 141|. Similarly, the 
phonon counting field ^(i) is equal to -l-(— )f on the pos- 
itive (negative) branch of the C-contour and stands for 
a virtual measurement of the energy transferred to the 
local phonon mode, i.e. the net variation of the number 
of vibrational quanta during the charge transfer process. 



The charge transfer mechanism in the molecular 
junction shown in FigU] is a coherent, time dependent 
process involving tunnelling of charges from the elec- 
trodes to the dot and energy exchange with the local 
phonon mode (emission and absorption of phonons). In 
the steady state, the electronic current (rate of charge 
transfer per unit of time) and the vibrational population 
are constant in average. However, both quantities might 
fluctuate in time. The most complete information about 
those fluctuations is encoded in the joined probability 
distribution P{q,n), defined as the probability that q 
charges are transferred through the junction and the 
number of phonons populating the mode has varied by 
an amount of n quanta during the measuring time io 
(this time is large enough to reach the stationary state). 
This quantity generalises the notion of full counting 
statistics (PCS) of transmitted charges which is recov- 
ered by tracing out the vibrational degrees of freedom 
Pei{q) = X^nGZ -f'C'?'"')- Similarly, one could obtain a 
phonon counting statistics by tracing out the electronic 
degrees of freedom Tph{n) = X^qgz ^(9' "■)■ ^"^ general, 
the joined distribution cannot be factorised due to the 
presence of e-ph correlations, i.e. P{q,n) ^ Pei{q)Tph{n). 

The aim of this article is to provide a way of com- 
puting this joined distribution by treating electronic 
and vibronic fluctuations on the same footing. By 
analogy with the electronic case [l|, y, [S^, we define 
the cumulant generating function (CGF) as the Fourier 
transform of the distribution P{q.,n), namely S'(x,C) = 
- InJEgez E„ez e''-i^+<^ P{q, n)}. This functional gen- 
erates the cumulants of the distribution P{q, n) by suc- 
cessive derivations of the CGF with respect to the elec- 
tronic (vibrational) counting field x (0- The connection 
of this quantity to the Hamiltonian of Eql2] is made by 
computing the following contour-ordered evolution oper- 
ator in Keldysh space 



C. Expansion of the CGF in the tunnel limit 



5(x,0 = -ln(Tcexp{ 



dtHT;x{t),i(t)]) ^ (3) 



In Eq|31 the tunnelling operator is obtained by perform- 
ing the following substitutions in the Ht Hamiltonian, 




I X=L,R 



FIG. 2: Linked cluster expansion of the CGF ^^^^^(x,^ to 
lowest order in V. Plain (daslied) line corresponds to the 
dot (electrode) Green function. The wavy line stands for the 
phonon Green function. 



Evaluating Eql3] for any value of the parameters 
(r,A^/ajo) is a formidable task. As a particular case 
of our general framework, we would like to address in 
detail the tunnel regime, for which F <C A^/a;o,wo and 
the charges are injected sequentially into the dot. Our 
calculation holds for off-resonant (low-transmission) 
situations, for which jwrfj ^ F and correlated transmis- 
sion events (electron bunching) may be neglected. In 
this regime, the resulting electronic FCS is anticipated 
to be Poissonian, but as far as we know, the explicit 
analytical evaluation of the FCS is still not reported in 
the literature. 

We thus perform a linked cluster expansion of the CGF 
to lowest order in the tunnelling rate F, or in an equiva- 
lent manner to the lowest order in the transmission coef- 
ficient (see Figl2]). Even if restricted to order F, this cal- 
culation holds for arbitrary values of the e-ph coupling 
strength g. The corresponding expression of the CGF 
provides in a single, compact formula, the most complete 
information about both electronic and vibrational fluc- 
tuations in the tunnel limit 
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+ W2) (4) 



where the counting field dependent tunneUing matrix 
element T^^^].^ = ri(j^)e±'("2-"i)x/4 takes into ac- 
count charge transfer processes from the left (right) elec- 
trode to the dot. Evaluation of Eq|4] involves the bare 
dot Green function gdd{t) = —^{Tcd{t)d^O))o (plain 
line in FigH]), the X = L,R lead Green function 
gxxit) = -i{T,^xit)^U0))o (dashed line in FigH]) 
and the ^-dependent phonon Green function V^{t) = 
(TcX(i)Xt(0))o (wavy line in FigE]). Exphcit expres- 
sions for those Green functions are given in Appendix 
A. 



III. RESULTS IN THE TUNNEL REGIME 
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A. Electronic FCS 



We derive from Eq|4] an analytical expression for the 
electronic CGF in the tunnel regime, namely SS^^ (x) = 
<55(i)(x,e = 0) 



ss. 



l^\x) = -to{rL- 



.R{e'^-l)+rR^L{e- 



1)} (5) 



Its Fourier transform provides a bidirectional Poissonian 
distribution for the electronic FCS Pei{q) with corre- 
sponding left to right (right to left) rates Tl^r (Tr-^l)- 
Those coefficients are evaluated within the scope of an 
approximation derived in Appendix B, which is by con- 
struction current conserving and consistent with the non- 
interacting limit g —i- 0. We obtain in the case of sym- 
metric coupling to the leads F^ = Tr = T/2 



FIG. 3: Current-voltage characteristics I{V) in the scope 
of the approximation derived in Appendix B. The parame- 
ters used for this plot are F = O.OOlojo, T — O.OSloq and 
fiL = —fJ'R = V/2. Black, blue and red curves correspond re- 
spectively to : A = 0.0; 0.3; l.Owo. Dashed, dashed-dotted and 
plain curves are obtained respectively with the dot position 
Oja = 0.00; 0.20; 0.75a;o. 



does not contain interaction-induced corrections to the 
Poissonian distribution. Such corrections were shown to 
arise from an avalanche (electron bunching) mechanism 
for the dynamics of the transferred electrons [l3| which is 
not accounted for in the lowest order expansion presented 
in EqlH Such terms will emerge from a non perturbative 
evaluation of the CGF ( "all order in F" dressing of the 
dot Green function), which is out of the scope of the 
present article [42 1. 



Tl^r - -{fLTRh + il-fR)TLe}{':^d) (6) 

Ffl^L = \{fRTLh + {l-fL)rRe}i<:^d) (7) 

In Eq|6] and [71 the electronic and hole rates for multi- 
phonon processes in the lead X = L, R are respectively 
defined as Fxe(w) = Fe^^ Y,n=o^fx('^ + "'^o) and 
Txhiio) = Te-a' j:tZ^{^ - fx{oJ - nwo)]. It is in- 
teresting to notice that the left to right and right to left 
rates of the Poissonian distribution exhibit a non ana- 
lytical behaviour in the e-ph coupling strength g at low 
temperature, i.e. a non perturbative calculation in g 
is necessary to derive correct results in the tunnel limit 
F — )■ 0. Compared to the existing results in the litera- 
ture 12l4l4|. the expression for the CGF derived in Eq[5| 



1. Electronic current 



The first cumulant < q > 



(1) = 



-(I) 



^rL,l corre- 



sponds to the mean current liV) that flows across the 
junction. We show on Fig[3| typical I{V) curves ob- 
tained by varying the dot position and the e-ph coupling 
strength. For the case aid = 0, the I{V) characteristics 
exhibit an inelastic threshold at T^ = i2LUQ corresponding 
to the activation of inelastic tunnelling, i. e. an electron 
on the dot may tunnel to the leads by emitting a phonon. 
Additional inelastic channels open when increasing A for 
voltages multiple of y = ±n2ujo, and correspond to the 
onset of multiple phonon emission. It is interesting to 
notice that those inelastic patterns are simply explained 
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FIG. 4: Noise-voltage characteristics S{V) in the scope of 
the approximation derived in Appendix B. The parameters 
used for this plot are F = O.OOlojo, T = O.OStiJo and ^l = 
— /ifl = V/2. Black, blue and red curves correspond respec- 
tively to : A = 0.0; 0.3; I.Owq. Dashed, dashed-dotted and 
plain curves are obtained respectively with the dot position 
(id = 0.00; 0.20; 0.75ajo. 



by Pauli principle that forbids multi-phonon transitions 
if the final channel of diffusion is already occupied. In the 
low-temperature limit F <^ T <C cjq, the shape and mag- 
nitude of the jurnps at V^ « ±n2u!o is the result of Franck- 
Condon factors [ll|,ll9| entering into Eql6]and[7l Strictly 
speaking however, our approximation breaks down in the 
limit ujd ^ (dashed curves of Figl3]) corresponding to 
resonant tunnelling. Only far from a resonant situation, 
namely when \ujd\ ^ F, is our low-F (and thus low- 
transmission) approximation valid (dashed-dotted and 
plain curves of FiglSj. As explained in Appendix B.l, 
the transmission factor of the non interacting junction is 
approximated in this limit by T{uj) w ■kT6{ll: — uJd), and 
the corresponding I{V) curves behave in the low-voltage 
region (0 < w^, F <C cjq) as 



<q> 



(l)c 



io- 



Te-3 



{fLi<^d) - fRi<^d)} 



(8) 



EqJH] is consistent with the Landauer-Biittiker formula 
331 , but with a renormaliscd rate of tunnelling F = Fe"^ 
43| . It fully coincides with it in the non- interacting case 

(g — >■ 0) and in the appropriate low-F limit stated above. 



2. Electronic current noise 



A similar behaviour (presence of inelastic thresholds) 
is observed for the case of the second cumulant (q^)c = 
F 



in 






that corresponds to the current-noise char- 
acteristics S{V) on FigUl In the case of low-voltages 
(0 < a)rf, F ^ Wo), the noise is given by 



<q'>i'^^to 



re-9 



{fL{l-fR) + fR{l-fL)}{':^d){9) 



In the appropriate low-F limit, EqlH] coincides with the 
scattering result for shot noise [6| when 5 — )► 0. 



B. Population dynamics of the phonon mode 

1. Derivation of the master equation for the phonons 

The computation of the phonon stationary distribu- 
tion, namely the probability Pph (n) of having n phonons 
populating the mode in steady state, is more difficult to 
achieve than the corresponding calculation of the elec- 
tronic FCS. The source of the difficulty originates from 
the low-F expansion of Eq|4] that docs not take into ac- 
count phonon emission and absorption on the same foot- 
ing, i.e. EqU] provides multi-phonon emission processes 
at this order in the F-expansion but lacks higher-orders 
absorption processes that are necessary to reach a steady 
state. 

We adopt in the following a self-consistent treatment 
that cures the problem by including in Eq|3] any excited 
state \m) of the local vibrational mode. More generally, 



n{m). 



we define TK[''{n) = X^ogz -^ (?' "") ^^ ^^® probabil- 
ity of transition from the initial vibrational state \m) 
to the final state \inii -t- n) during the measuring time 
to, after tracing out the electronic degrees of freedom. 
This quantity is computed from the generalised vibra- 
tional CGF, namely 5S^^l^{i) = 5sln\x = 0,0, that 
includes the contribution of vibronic Green functions 
V^^\t) = (rcX(i)Xt(O))™ averaged other the excited 
phonon state (to quanta in the phonon mode). A de- 
tailed derivation of the generalised vibronic Green func- 
tions and phonon CGF is proposed in Appendix A and 
C respectively. 

In the absence of any external damping mechanism for 
the phonons (non equilibrated phonons), the dynamics of 
the phonon population is self-determined for each time 
interval t^ by the electronic tunnelling mechanism as 



Pph {to ; 



-I- 00 

m=0 



{n~m)PphiO;m) (10) 



In Eq llOl Pph{ta;n) is the probability of having n 
phonons populating the local vibrational mode at time 
to. The stationary distribution of the phonons is ob- 
tained in the long time limit as the fixed point of Eq llOl 
In the tunnel limit F ^- 0, the transition rate Ip™ (n) is 
expanded in leading order of F as 



T^in) 



Sn,o + tori';:\n) + o{tlT') 



n(m) 
ph 



in) 



' ph 



iOe 



-in^ 



(11) 

(12) 



^(m) 



In Eq llll rL (n) is a transition rate per unit of time cor- 
responding to the multi-phonon process \m) — >■ \m -\- n). 
It is formally related to the Fourier transform of the 

It is interesting to 



phonon CGF S'J"^(n) (see EqlE]) 
notice that the matrix \rph\n,m — T. 



m) might 
be evaluated analytically (see Appendix C) and has the 



(m)i 
ph 



property of conserving the normalisation of the phonon 
distribution, i.e. for each index m we have the follow- 
ing relation amongst matrix elements '}2,n=o ^ph:n,m = 0. 
Eqll2land ll3l arc an important result of this paper. They 
are the constitutive relations that connect the formulation 
of transport based on perturbation theory in Keldysh space 
(usually used in the coherent transport regime) to the one 
based on the master equation in the tunnel limit. Includ- 
ing Eqllll into EqllO| we derive from our formalism the 
standard quantum master equation for the dynamics of 
the phonon population 
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m— 0;T7i7^n m—O-^m^n 



(13) 
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2. Results for the phonon stationary distribution 
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FIG. 5: Stationary distribution of the phonon population 
Pph{n) for a symmetrically biased molecular junction with 
parameters F = O.OOlwo, T — 0.03a;o and V = 2.5a;o. Blue, 
orange and red curves are obtained for : A = 0.3; 1.0; 2.0 wo 
respectively. Plain, dashed and dashed-dotted curves corre- 
spond respectively to the dot position CJd ~ 0.00; 0.20; 0.75a;o. 
The black curve corresponds to 'F = O.Sojq, ij^d ~ 0.00 and 
A = O.SojQ. 

As an example, we compute the phonon stationary dis- 
tribution Pph{n) = \irat^+oo Pph(t]n) obtained as the 
zero eigenvector of the [F^/j] matrix (see Ea fT3| . For 
voltages below the inelastic threshold {V = O.Swq for 
the black curve of FiglS]), phonon emission is forbidden 
by Pauli principle and the phonon mode is not populated 
whatever the strength of e-ph coupling. In this case, the 
distribution Pph{n) is a peak at n = given by the Bose 



equilibrium distribution. For voltages above the inelas- 
tic threshold, a pumping mechanism appears : the on- 
set of phonon emission strongly drives the population of 
the vibrational mode out of equilibrium. For weak e-ph 
coupling (A ~ O.Swq) and V = 2.5ajo (see blue curves 
on FiglS)), the distribution Pph{n) is long tailed, whereas 
for increasing A = 1.0 — 2.0loo (see orange-red curves on 
FigEI, it gets closer to the equilibrium distribution (the 
tail is shorten). The corresponding non- monotonous be- 
haviour of Pph{n) with A is resulting from the compe- 
tition between the emission-absorption mechanism given 
by EqUni and the selection rules imposed by the tran- 
sition matrix [Fp/j] (when A increases, desexcitation of 
highly-excited states toward the ground state are more 
likely to happen[19[). 



CONCLUSION-OPEN QUESTIONS 

In this work, we have developed a theoretical frame- 
work that enables to treat on the same footing the 
concept of counting statistics of transmitted (fermionic) 
charges commonly used in the field of quantum trans- 
port and the concept of counting statistics of phonons 
(bosons). We illustrated this framework by studying 
in detail the transport properties of an interacting 
molecular junction, from both point of view of counting 
statistics of transmitted electrons and excited phonons. 
We derived non-perturbative results in the e-ph coupling 
strength which are valid in the lowest order of the 
tunnelling rate to the leads F (tunnel regime). Our 
description provides a general framework that unifies 
previous studies based on rate equations (in the tunnel 
regime) to the one based on non-equilibrium Green 



functions techniques (in the coherent regime). 

A natural extension and open question arising from 
the present work is to investigate the role of higher- 
order terms in the F-expansion of EqU) In particular, 
co-tunnelling processes [SJ, |35| are expected to be im- 
portant in order to describe off-resonant transport into 
the deep Franck-Condon blockade [9| or to derive results 
valid in the limit of large transmission (resonant tun- 
nelling) [36|, |37| . Following the recent work of Maier et 
al |38| . such an issue could be investigated by using a 
self-consistent scheme to compute the FCS in the strong 
coupling regime. More generally, the nature of trans- 
port properties and mechanisms in the polaron crossover 
regime is an important, unsolved problem [39[ that de- 
serves future investigations following the lines drawn in 
the present article. 



9dd(.^) = X {^) (15) 



9dd i^) = «r V (cj) 



(16) 



-+( ^ .^ e-^^/^(l-/L)+e»^/^(l-fe) , ,^.. 

5d/('^) = -*r (wjji?) 

where the Keldysh determinant is given by 



1\^{LU)^{L0-Ujdf+T^+ (18) 

r'Ke^'^ - 1)./l(1 - Ir) + (e-^^ - l)/fl(l - /l)}(w) 

The mean population of the dot < n^ > is obtained af- 
ter integration of the non-diagonal component of the dot 
Green function g~^^{^) expressed at zero counting field 
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< nd >= 



duj 



T 



2Tr{uj-LJdf+r^ 



{hiio) + fR{u)} (19) 



In the atomic limit (F ^> 0), far from resonance {\ujd\ ^ 
F), the Lorentzian-integrand of EgllDlis approximated by 
a delta function and we obtain for the voltage dependent 
mean population of the dot 



lim < Ud >= 7r{/L(wd) + /i?Xwd)} (20) 

r->o 2 

Similarly, the atomic limit of the non-interacting Green 
function is obtained as 
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Appendix A : free Green functions 



1. Electronic Green functions 



gZH - «{ 



I- <nd> 



< nd > 



(jj — (jjd + onrj uj — ujd - 
9ddi^) = *27r <nd> 5{uj - Ud) 
9ddi'^) = -i27r(l- < Ud >)5{uj - ujd) 



-} (21) 
airj 

(22) 

(23) 



In this part, we compute the non-interacting electronic 
Green functions of the dot and of the X ~ L,R leads, 
namely gdd{^) and gxx{^)- 

We suppose that the leads are maintained in equilibrium 
and characterised by a flat density of states po and a 
chemical potential /-iL(fl) = +{~)V/2. The correspond- 
ing free Green functions are given by 



gxxi^) 



W 



2/x - 1 2/x 

-2(1 -/x) 2/x-l 



(u) (14) 



where /x(w) is the Fermi distribution of the lead X = 
L, R and W = I/ttpq its bandwidth. 
The free (non interacting) dot Green function is deter- 
mined in the case of symmetric contacts to the leads 
Tl=Tr = F/2 as 



2. Vibronic Green functions 

The free vibronic Green function T'(^=o{t) ~ 
{TcX{t)X''{0))Q is evaluated by performing the aver- 
age other the phonon ground state {m — quanta in 
the vibrational mode). Using the Glauber equality |3| 
e'^'^^ — e'^e^e^al'^^^l which is valid for any operators 
(C, D) that commute with their commutator, we obtain 
similarly to Ref. [40| 



n"=o(^) 



r?-o^Lo) 



1 



+00 2ti 

^-^ n\ oj — anujQ + iri 

n— 

1 



} 



— -2n6(uj + anujo) 



Lo + anojQ — ir] 

+°° 2n 

n 



(24) 



(25) 



n=0 



More generally, we define the ^-dependent free vibronic where the correlators C"^(t) = {m\Xa^{t)XljO)\m) 
Green function as V^"'\t) = (rcX(i)Xt(0))„, where the and C"^.„(i) = {m\xlAO)Xfi^{t)\m) are given in Fourier 



mean value is performed other an excited phonon state 
(m ^ quanta in the mode). We obtain in general 



representation by 






0{t)C:^.n{t)+0{-t)CZ%,r.{t) (26) 

CZ+{t) (27) 



-tQ/3 



Cri.,,ni^)^2ne' 



2^ 2^ 2^ Ti, i^„i ""P2""P3 

Pl=0p2=0p3=0 



PlW- 



c-Cp^r^-ire- 



-(Pl+P2-P3 



^^S[uj =F [pi +P2 -ps) 



1^0 



(28) 



It is interesting to notice that only non-diagonal compo- 
nents of the phonon propagator are explicitly dependent 
on the phonon counting field. If the phonon states are 
restricted to the to = ground state (like in the low-F 
result of Eq|4]) , the only possible processes available cor- 
respond to multiple phonon emission. Only when con- 
sidering an excited phonon state {m > 1) are absorption 
processes allowed. 



Appendix B : rates of the Poisson distribution 

1. The non interacting case 

Using our formalism and the expression of the dot 
Green functions (see Eg lTMTT)) . the electronic CGF of 
the resonant level in the non-interacting case {g = 0) 
might be determined in all orders of F. This provides 
the Levitov-Lesovik result [l|, [2, |32| for the free electronic 
CGF S'g; (x), and the associated binomial FCS of trans- 
mitted electrons 



the Poissonian result S^^ (x) 



to{ri°!,^(e^'^ 



sTix) 



-to 



/^,„{:.., 



(29) 



(e^^ - 1)/l(1 - /fl) + (e-*'^ - l)/i,(l - h)] H} 

Eql^ni is associated with resonant tunnelling of electrons 
and holes through the dot, characterised by a transmis- 
sion factor 



T{u) 



[U - LUdY + F2 



(30) 



The limit of low-F and off-resonant situation \ijjd\ 3> 
F corresponds to an approximate transmission factor 
T{oj) K, ■kT5{lj — ujd)- The corresponding expansion of 
the electronic CGF in the lowest order of F provides 



-.(0) 



{e 



«x 



1)} with rates expressed as 



f2,H = \fL{u^a)[l-fE{u:a)] 



1) + 

(31) 

(32) 



In the following, the results for the interacting case have 
to be understood as well in the appropriate low-F and 
low-transmission limit stated above. 



2. Interacting case : choice of the gauge 

When evaluating the lowest order in the F-expansion 
of the electronic CGF, one is not ensured (if no self- 
consistency is achieved) that current conservation is ful- 
filled, namely that the obtained cumulants versus voltage 
curves have a well defined symmetry under the transfor- 
mation V ^ —V, nor that the obtained CGF coincides 
with the non- interacting result of Eg I5l1 and Eg 15^ when 

We use this constraint of both charge conservation and 
recovering the non-interacting Levitov-Lesovik result to 
select the proper gauge in implementing the electronic 
counting field in the Hamiltonian. Clearly, the choice of 
a gauge that incorporates the electronic counting field in 
the left electrode only tLd — >■ iidS"*'^'*'/^ breaks the sym- 
metry between left and right electrode when expanding 
the CGF to the lowest order in F, i.e. this choice of gauge 
is not current conserving. The symmetry between both 
electrodes (and hence current conservation) is restored by 
introducing the electronic counting field symmetrically in 
both leads, namely iL(_R)d — >■ iL(fl)de^*-^^*^^''. In the fol- 
lowing, we make this choice of a symmetric gauge that is 
by construction current conserving. 
We then evaluate the electronic CGF from EqHl using 



Ea l24ll25] for the free phonon propagator 



+ °° „2n 



n^ — n ■ ^ L -J 



OL — OL / \ 

9dd ('^i) 



r^e 



-*"x/2 



9ll +ri?e 



Jax/2 



ffiii? 



(HJI 



anojo 



)} (33) 



r 



In Eg 1551 the free dot Green function gdd{^) has to be 
evaluated. However, by implementing naively the bare 
dot Green function given in the atomic limit by Ea l21f 
[23l the obtained CGF although current conserving does 
not reproduce the limiting case of the Levitov-Lesovik 
formula when g -> 0. This is related to the fact that 
the non equilibrium state of the dot level subsystem is 
ill defined in the atomic limit T — > 0, i.e. one has to 
artificially include the presence of electrodes maintained 



under a constant voltage bias as an external boundary 
term (for instance, as a voltage dependent population of 
the dot in the expression of the bare dot Green function 
in Eq l2Tll23| . We have found the following procedure to 
overcome this difficulty and compute safely the electronic 
CGF in such a way that the non-interacting limit for the 
rates as given by EqEUEH is recovered. We first use the 
expression of the free dot Green function in all orders of 
r as written in Eq fTMTTl and obtain for the CGF 



J 



fRALh + {l-fL)AR,]}{Lj) 



(34) 



IlArh + (1 - .fB.)ALc 



We finally evaluate the Keldysh determinant A^(w) at 
zero electronic counting field x = and go to the limit 
r ^ 0. We obtain EqUKZlfor the electronic CGF and its 
corresponding rates. 



sponding multi-phonon transition rates per unit of time 

^p"i ("")■ ^^ -^^^^ write EqHlin terms of the generalised 
phonon Green functions defined in Appendix A 



Appendix C : transition rates for the phonons 



In this Appendix, we derive analytical expressions for 
the generalised phonon CGF ^5*^,^ „j(^) and for the corre- 



+00 m 2p2 2Cni+no~l 



Pl=0p2=0p3=0 



PiW "' "' 



fjmfj2p2 (—\\P3 



J{Pl+P2-P3)i 



(35) 



Yl / -Kr{9dd"i^i)[^L9LL" + ^R9rT {uJi - a{pi + P2 - P3)^q)} 

a=±'' 



r 



Putting the evaluation of the leads and of the dot Green Fourier transforming Eql35] the following expression for 
functions as given by EalTiland Eq l^TII^ wc obtain after the transition rates 
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^(m) 



rrw 



+ 00 m 2p2 

-^^ E E E 



2(P1+P2) 



Pl=0p2=0p3=0 

<nd> [Tl{1 - II) + Tr{1 - fR)]{Lo^) + 



Pi!p2! 



I 



^n.pi- 



-P2-P3 



^n,0 



(36) 



where a;± = (Dd ± (pi + P2 — P3)wo- The constraints on 
summations due to the last term in EqlHSl imply that the 
transition rates are conserving the phonon probabihty, 



namely that X^n^ ^'. 



ph^n^ni 



0. 
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